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ABSTRACT 

NRL  has  had  several  artificial  satellites  placed 
in  orbit  for  various  experimental  purjxjses.  By  ob¬ 
serving  the  changes  in  orbital  elements,  it  lias  been 
possible  to  calculate  the  density  at  the  altitude  range 
of  950  to  1160  km  for  the  past  3  years.  Long-  and 
short-term  effects  are  noted  and  discussed.  A  cycli¬ 
cal  variation  of  27  days  has  been  found.  It  is  appar¬ 
ently  due  to  the  variation  in  solar  activity  with  the 
sun’s  rotation.  The  centers  of  activity  appear  to  be 
stationary  on  the  sun's  surface  for  months  at  a  time. 
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UPPER-ATMOSPHERE-DENSITY  DETERMINATION* 
USING  NRL  ARTIFICIAL  SATELLITES 


NRL  has  had  several  satellites  placed  in  orbit  for  various  experimental  purposes. 

As  a  part  of  the  development  of  the  Space  Surveillance  System  (1),  several  smooth 
spheres  have  been  launched.  Observations  of  these  objects  were  used  to  determine  ti:e 
system  coverage.  Two  identically  appearing  spheres,  having  a  10:1  ratio  in  masses, 
formed  the  basis  for  the  Dual  CalSphere  experiment  (2).  The  changes  in  semimajor 
orbital  axis  of  these  two  sphc:  es  were  recorded  to  calculate  fie  relative  drag  effects  of 
the  two  objects.  The  rate  of  change  of  the  semimajor  axis  was  seen  to  vary,  even  in  the 
few  weeks'  time.  It  appeared  that  extended  observations  of  the  lighter  calsphere  should 
provide  upper-atmosphere-density  data  for  as  long  a  period  as  one  wished  to  observe. 

With  this  in  mind  along  with  the  need  for  a  large  target  satellite  for  Space  Surveil¬ 
lance  System  calibration,  NRL  built  and  arranged  for  the  launch  of  two  dodecapole  re¬ 
flectors.  Each  object  is  a  26-cm  sphere  with  12  hollow  rods  erected  normal  to  the  faces 
of  an  equivalent  regular  dodecahedron.  Each  rod  is  1.27  cm  in  diameter  and  7.62  m  long. 
A  photograph  of  the  object  with  rods  not  extended  is  included  as  Fig.  1.  These  satellites 
serve  several  purposes.  To  the  217-mHz,  bistatic  cw  radar  of  the  Space  Surveillance 


Fig.  1  -  Dodecapole  satellite  (prelaunch) 
with  the  12  rods  in  stowed  position 
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System,  they  are  large  targets.  For  the  same  reasons,  they  are  useful  as  passive  com¬ 
munication  reflectors.  Finally,  the  mass-to-area  ratio  and  the  nearly  isotropic  projected 
area  make  them  useful  objects  for  determining  drag  due  to  the  atmosphere  at  the  1000-km 
altitude. 

These  objects  and  the  facilities  of  the  Space  Surveillance  System  make  data  available 
for  a  region  of  the  upper  atmosphere  thzt  is  otherwise  poorly  covered.  Little  photo¬ 
graphic  work  is  done  on  objects  at  this  height. 

As  a  part  of  its  regular  task,  the  Space  Surveillance  System  prepares  five-card 
orbital  elements  on  over  a  thousand  orbiting  objects.  Current  observations  are  used  in  a 
differential  correction  program  to  update  the  elements.  The  semimajor  axis  is  calcu¬ 
lated  to  the  eighth  decimal  place  in  terms  of  the  earth’s  radius.  In  the  appendix,  the 
change  in  semimajor  axis  is  related  to  the  density  of  the  medium  in  which  the  satellite 
moves. 


Table  1 

Excerpt  from  Log  of  Significant  Data  on  Object  900 
(Calsphere)  and  Object  1310  (Dodecapole) 


Date 

(Modified  Julian  Days) 

Date 

(Calendar) 

Semimajor  Axis 
(Earth  Radii) 

Object  900 

Object  1310  | 

39959 

Apr.  13,  1968 

1.16717726 

1.14241144  | 

39966 

Apr.  20,  1968 

1.16717258 

i.  14238240  | 

39973  j 

j  Apr.  27,  1968 

1.16716880 

1.14235698  | 

The  log  is  continuous  from  January  1966  to  date.  The  elements,  with  epoch  at  the 
last  south-to-north  equator  crossing  on  Friday,  are  received  at  NRL  each  Monday  morn¬ 
ing.  The  density  values  for  each  week  are  plotted  in  Fig.  2  and  connected  by  segments  of 
straight  lines.  The  smooth  curve  joins  points  representing  the  mean  value  of  each  week 
and  the  two  before  and  the  two  after.  These  5-week  running  times  tend  to  obscure  short  - 
term  variations,  1  month  or  less,  but  point  up  the  longer  period  trends.  The  lower  plot 
is  for  satellite  1964-63C,  the  lighter  sphere  used  in  the  Dual  Calsphere  experiment.  The 
upper  plot  is  for  1965-16G,  one  of  the  dodecapoles. 

The  weekly  changes  in  semimajor  axis  are  shown  in  Fig.  3  for  three  different  satel¬ 
lites.  Both  dodecapoles  and  the  light  calsphere  are  shown  here  with  a  logarithmic  scale. 
The  purpose  is  to  show  the  close  time  correlation  of  the  accelerations  of  the  three  ob¬ 
jects.  The  plot  shows  that  the  ratios  of  the  maximum -to-minimum  acceleration  of  the 
sphere  and  dodecapole  at  nearly  the  same  height  are  approximately  equal.  The  dodeca¬ 
pole  at  lower  altitude  shows  a  somewhat  greater  ratio. 

Looking  at  Fig.  2,  it  is  apparent  that  the  air  density,  as  derived  from  acceleration, 
has  shown  marked  changes  during  the  last  2-1/2  years.  Added  to  the  general  increase 
are  at  least  two  apparently  periodic  functions.  The  weekly  data  points  show*  13  maximum 
values  during  1967.  The  mean  period  is  approximately  27  days.  This  interval  is  close 
to  the  sun's  synodic  rotation  period  and  also  the  moon’s  period.  The  smooth  curves  show 
maximum  values  some  10  months  apart.  Note  that  the  short-term  variation  has  maxi¬ 
mum  amplitude  at  the  10-month  maximum  times.  This  observation  tends  to  eliminate 
the  influence  of  the  moon,  which  Lr  significant,  would  be  constant  except  for  a  small  in¬ 
teraction  with  the  sun.  The  27-day  variation,  then,  must  be  attributed  to  solar  activity, 
varying  with  rotation. 
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Fig.  3  -  C  tangos  in  so  mi  major  axes  of  three  satellites  wrsss  ticso 


The  variation  is  consistent  with  the  effect  of  a  source  of  solar  activity  —  stationary 
on  the  sun’s  surface  and  varying  in  intensity  over  a  period  of  months.  The  activity  ap¬ 
pears  to  persist  through  the  changing  pattern  of  sun  spots. 

At  this  time,  November  1968,  accelerations  are  high,  indicating  the  third  such  period 
of  maximum  activity  since  the  start  of  the  study.  The  apparent  10 -month  period  is  sup¬ 
ported  by  these  three  maxima,  but  further  observation  -ill  be  necessary  before  any  firm 
inference  can  be  drawn. 

All  calculations  made  to  this  point  are  based  cm  the  assumption  that  each  weekly 
decrement  in  the  semimajor  axis  reflects  a  constant  or  slowly  changing  drag  during  that 
period.  Some  variations  must  be  considered  to  see  if  the  assumption  is  justified.  The 
effective  density  of  the  upper  atmosphere  is  a  function  of  the  temperature  resulting  from 
solar  excitation-  Durieg  this  study  each  satellite  has  ueen  subject  to  various  geometrical 
arrangements  with  respeci  to  the  earth  and  sun.  At  times,  each  object  was  in  sunlight 
throughout  the  orbit.  At  other  times  they  had  varying  times  of  immersion  in  the  earth's 
shadow.  At  times  the  plane  of  the  orbit  included  the  diurnal  solar  bulge.  Since  the  ob¬ 
servations  cover  more  than  2  years,  seasonal  variation  in  sun-to -earth  distance  must  be 
considered. 

If  the  earth -sun  line  is  nearly  perpendicular  to  the  plane  of  the  orbit,  the  satellite 
will  move  through  a  medium  that  is  uniformly  excited  by  the  sun.  At  the  altitude  of  these 
satellites  the  surface  temperature  variations  have  little  effect  on  the  atmospheric  tem- 
peratur ;.  As  the  plane  of  the  orbit  moves  to  include  the  sun’c  position,  the  satellite  en¬ 
counters  regions  of  different  excitation  including  some  fuli  shadow.  If  the  integrated 
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effect  of  the  drag  around  the  stole  orbit  varies  as  a  function  of  the  apparent  solar  posi¬ 
tion,  then  a  periodic  effect  should  be  noted.  Since  the  rotation  of  the  plane  of  a  near 
polar  orbit  (satellite  900)  is  only  a  few  degrees  per  year,  the  effective  period  of  any  var¬ 
iation  tfete  to  solar  position  should  be  nearly  1  year. 

Due  to  the  economic  problem  of  getting  satellites  launched  as  cheaply  as  possible, 
XRL  has  taken  advantage  of  excess  lift  capability  of  certain  primary  payload  operations, 
even  though  the  orbital  parameters  might  not  be  exactly  those  desired,  rbe  dodccapole, 
object  1310,  ns  placed  in  an  orbit  sitb  a  nominal  70-degree  inclination.  In  this  case, 
the  result  is  good.  The  orbital  plane  revolves  through  360  degrees  in  approximately  4 
months  as  compared  sitb  the  12-mocih  period  of  the  calspfeere  (900).  The  close  agree¬ 
ment  of  the  accelerations  of  the  tso  bodies  shore  that  the  effect  due  to  solar  position  Is 
small  cr  negligible  compared  with  that  due  to  charges  in  solar  radiation. 

Artificial  satellites  of  known  size  2nd  mass,  isotropic  drag,  and  circular  orbits  used 
in  conjunction  with  current  orbital  element  determinations  can  provide  good  measure¬ 
ment  of  the  density  of  the  outer  atmosphere.  The  objects  used  in  this  study  are  well 
fitted  for  the  altitude  near  1090  km.  The  accelerations  are  high  enough  to  measure 
easily  and  low  enough  that  data  may  be  obtained  over  an  extended  period  of  time. 

At  a  lower  aititude  these  objects  would  decay  so  rapidly  that  they  could  show  only 
short-term  effects.  This  points  up  the  need  of  some  future  launches  of  spherical  bodies 
with  greater  mass-to-area  ratio  for  lower  altitude  studies. 

An  object  with  higher  drag  could  be  used  with  a  daily  clement  updating  pc  ogram  to 
make  a  more  detailed  study  of  short-term  variations  that  is  possible  with  the  present 
weekly  reports. 

This  XRL  study  will  be  continued  through  the  coming  period  of  higher  solar  activity 
and  as  lotsg  thereafter  as  the  data  arc  significant. 
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Appendix 


ATMOSPHERIC  DENSITY  AT  600  SAUT  MI  AS  DETERMINED 
BY  THE  CHANGE  OF  ORBITAL  OF  SATELLITE  900, 

ONE  OF  NHL  CALSPKERES 

The  density  of  I  be  upper  atmosphere  has  been  calculated  in  several  studies  froen 
observations  of  satellites.  There  is  a  direct  relationship  between  the  decay  off  a  satellite 


The  density  of  the  upper  atmosphere  has  been  calculated  in  several  studies  froen 
observations  of  satellites.  There  is  a  direct  relationship  between  the  decay  off  a  satellite 
orbit  and  the  air  density.  In  renera],  the  rate  of  change  of  the  senisajor  axis  of  an  ob¬ 
ject  in  orbit  about  the  earth  is  proportional  to  the  density  of  the  atmosphere  through 
which  it  movesf 
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where  <-  =  eccentricity.  Before  defining  the  other  elements  in  the  equation,  we  simplify 
it  by  setting  e  -  e?  therefore. 


,  J^SsLv 


k  -  **.  r 


where 


f  s  Telocity  of  the  satellite  with  respect  to  the  atmosphere, 
i't  =  velocity  of  the  satellite  with  respect  to  inertial  space, 
y,  *  velocity  of  the  atmosphere  with  respect  to  inertial  space,  and 
fi  -  angle  between  r,  and  the  plane  of  the  orbit. 

Assume  p  =  9Q  degrees  fore-0  and  therefore. 


rtrsr, 

-2-rr_ 

E2TT  * 


where 


f  -  atmospheric  density, 
CB  =  ballistic  coefficient. 
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*  -  projected  area, 

»  *  Mgshr  Telocity,  aad 
»  =  mss  of  satellite. 


All  of  these  parameters  are  easily  dtSstd  except  for  the  ballistic  coeOieiess  C&.  Is 
the  early  years  of  cpper  atgso^fetre  research,  the  raise  2.0  «as  ass— — ed.  Accord; sg  eo 
a  reeest  steady*  the  raise  for  as  cedde-ecated  sphere  is  the  a&ilafe  n^e  «bere  helisas 
is  the  principal  composes:  of  the  atmosphere  lies  ia  the  ra^ge  2.4  to  2.6.  For  a  polar 
octal  *"  5‘,;  therefore. 


f4C&ri 


where 


s'-  -  r »  j  . 


ia  which  =  curcsiar  orfcital  velocity  at  the  surface  of  the  earth  =  7.205  Iss/sec  aad 

l  g  -  a*  *  S»?  2-  . 
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The  ranges  sszshers  ia  the  shore  ccerersios  pre  a  fotssk  that  accepts  the  qwas- 
titles  i  aaf  »  ia  earth  radii  and  earth  nSi  per  week.  respectively.  Is  this  steady,  the 
raise  of  she  sessstjtx  axis  is  acted  each  week,  aad  the  dlftrtscc  betxtce  saccissrre 
readings  is  taxer  as  the  seas  rate  for  the  week. 
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